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ABSTRACT 

We combine a series of high-resolution simulations with semi-analytic galaxy formation mod- 
els to follow the evolution of a system resembling the Milky Way and its satellites. The semi- 
analytic model is based on that developed for the Millennium Simulation, and successfully 
reproduces the properties of galaxies on large scales, as well as those of the Milky Way. In 
this model, we are able to reproduce the luminosity function of the satellites around the Milky 
Way by preventing cooling in haloes with Fvir < 16.7 km (i.e. the atomic hydrogen cool- 
ing limit) and including the impact of the reionization of the Universe. The physical properties 
of our model satellites (e.g. mean metallicities, ages, half-light radii and mass-to-light ratios) 
are in good agreement with the latest observational measurements. We do not find a strong 
dependence upon the particular implementation of supernova feedback, but a scheme which 
is more efficient in galaxies embedded in smaller haloes, i.e. shallower potential wells, gives 
better agreement with the properties of the ultra-faint satellites. Our model predicts that the 
brightest satellites are associated with the most massive subhaloes, are accreted later (z < 1), 
and have extended star formation histories, with only 1 per cent of their stars made by the 
end of the reionization. On the other hand, the fainter satellites tend to be accreted early, are 
dominated by stars with age > 10 Gyr, and a few of them formed most of their stars before the 
reionization was complete. Objects with luminosities comparable to those of the classical MW 
satellites are associated with dark matter subhaloes with a peak circular velocity > 10 km s~^, 
in agreement with the latest constraints. 

Key words: galaxies: dwarf - galaxies: formation - Local Group - cosmology: theory - dark 
matter 



1 INTRODUCTION 

The satellites of the Milky Way (MW) are powerful touchstones 
for galaxy formation and evolution theories. Their proximity fa- 
cilitates detailed observations and characterisation of their prop- 
erties and hence constrains 'near-field' cosmology. In addition, 
their shallow potential wells make them more sens itive to astro- 
physical p rocesses such as supernova (SN) feedback l lLarsonI 19741 ; 
iDekel & Silk 1986) or to th e presence of a photoionization back- 
ground jBabul & Reeslll992h . 

Deep images have allowed the construction of colour- 
magnitude diagrams (CMD) of the MW satellites, from which the 
star formation histories have been deduced. These studies indi- 
cate that there is a large variety in the star form ation histories of 
these galaxies l lMateol 19981 : iDolphin et al.l2005h . The two gas-rich 
dwarf irregular (dirr) Magellanic Clouds show on-going star for- 
mation while the other dwarf spheroidal galaxies (dSphs) are gas- 
deficient and show generally little evidence for recent star forma- 
tion. Modem studies have revealed that all satellites contain an old 
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stellar population (> 10 Gyr), which likely keeps the imprints of 
the young Universe. 

In recent years, the number of known satellites around the 
MW has doubled, thanks to the discovery of very low surface 
brightnes s dwarf galax i es in the Sloan Digita l Sky S urvey (SDSS) 
Willman"e tai]|2005al l bl:lBelokurov et ai]|20 06. 2007V zucker et al ' 



2006a b; Irwin et all [20071 ; IWalsh et alj 12007. : .Belokurov et al 



2008, 200§). Since the sky coverage of SDSS DR5 is about 1/5 
of the full sky a nd the surface brigh tness limit is about /i ~ 
30 mag arcsec' -2 i IKodosov et ai]|2008l) . many satellites are likely 
yet to be discovere d in the next generations of surveys. E.g. 
iTollerud et al. 1 20081) have used sub-samples of the Via Lactea I 
( iDiemand et al. |2007|) subhaloes to conclude that the total num- 
ber of MW satellites within 400 kpc should be between ~300 and 
~600 and dominated by satellites fainter than Mv = —5. 

The new SDSS satellites have lower surface brightness > 
27 mag arcsec"^) compared to the classical satellites, but similar 
physical sizes. They have comparable luminositie s to some Galac- 
tic gl obular clusters, but are significantly bigger jBelokurov et al.l 
l2007h . The nature of these newly discovered satellites is still un- 
clear. They could be the prol ongation towards fainter luminosities 
of the classical MW satellites l lKirbv et alj2008h , tidal features (e.g. 
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Hercules dSph. lColeman et al .l2007l : ISand et al.l2009l) , or represent 
a completely new class of objects. 

Kinematic modelling based on line-of-sight velocity disper- 
sions, have demonstrate d that the classical MW satellites are dom- 
inated by dark matter jMateo et all 1 19931 : lMateo|[T993) . Recent 
studies have shown that, under the assumption of virial equilib- 
rium, the ultra-faint satellites have mass-to-light ratios as high as 
~ 100 — 1000, implying that th ese objects are the most dark matter 
domi nated systems known (e.g. iMurioz et alj200d : ISimon & Gehal 
l2007h . The constraint on the total mass inferred by velocity dis- 
persions is still uncertain because of the mass-velocity anisotropy 
degeneracy and of the small number of tracers employed in these 
studies. Recent analyses suggest that the MW satellites (including 
the newly discovered SDSS satellites) have a common mass scale 
when considering their in nermost regions within 600 or 300 pc 
dStrigari et all2007l , l2008l) . 

The cold dark matter (CDM) hierarchical paradigm suc- 
cessfully explains the large scale structures of the Universe 
jSpergel et al]l2007l) . Semi-analytic (hereafter SA) galaxy forma- 
tion models coupled with merger trees extracted from A^-body 
simulations, represent a useful technique to diagnose the complex 
physics involved in galaxy formation, with modest computational 
costs. In recent years, SA models have been proved to success- 
fully reproduce a number of observational measurements (e.g. spa- 
tial and colour-magnitude distributions) for galaxies seen in the lo- 
cal U niverse and at higher redshift (for a recent review, see lBaughl 
l200d) . In spite of the encouraging progress on the large scale, how- 
ever, CDM still faces severe challenges on the galaxy-scale and 
below. An example is the 'missing satellites problem' : namely the 
substructures resolved in a galaxy-size DM hal o significantly out- 
number the satelli tes observed around the MW jKlvpin et al.ll99^ ; 
iMooreetal .11 19991) . A number of studies have suggested that astro- 
physical processes such as the presence of a phot oionization back- 
ground m ight reconcile this discrepancy (e.g. Kauffmann et al.l 
Il993l : [Bul lock et al. 2000; Benson et al. 2002; Somerville 20oJx 
without invoking modifications on the nature of the DM parti- 
cles to reduce the power on small scales of the power spectrum 
jKamionkowski & Liddlell2000l ; IZentner & Bullocldl2003l) . 

Several groups have recently attempted to model the proper- 
ties of the M W and its satellites in a (semi-)cos mological setting. 
For example. lKravtsov. Gnedin & Klvpiiil | |2004|) have analysed the 
dynamical evolution of substructures in high-resolution A*'-body 
simulations of MW-like haloes and suggested that all the lumi- 
nous dwarf spheroidals in the Local Group are descendants of the 
relatively massive (~ 10® M©) high-redshift haloes, which were 
not significantly affected by the extragalactic ultraviolet radiation. 
iFontetalJliooa) have successfully reproduced the observed chem- 
ical abundance pattern of the MW stellar halo by combining mass 
accretion histories of galaxy-size haloes with a chemical evolution 
model for individual satellites. It should be noted that in these stud- 
ies, the phenomenological recipes adopted for star formation and 
feedback have been tuned to reproduce some of the properties of 
the s atellites in the Local Group. 

iBenson et alj | |2002|) have used a SA model which success- 
fully reproduces the present-day field galaxy luminosity function 
to study the properties of dwarf satellite galaxies known at the 
time. Their model calculates the influence of reionization self- 
consistently, based on the production of ionizing photons from stars 
and quasars, and the reheating of the intergalactic medium. This 
same model reproduces quite nicely the luminosity and size dis- 
tributions, gas content and metallicity of the classical satellites of 
the MW and M31. These authors have suggested that the surviv- 



ing satellites are those which formed while the Universe was still 
neutral. This study was carried out before the boost at the faint 
end of the satellite luminosity function. Although they did predict 
a large number of faint satellites below Mv ~ —5, after extrapo- 
lating their pr ediction on the luminosity-size space to the faint end, 
iBenson et alt s faint satellites tend to be too small at a given abso- 
lute m agnitude compared to the ultra-faint satellites discovered in 
SDSS i IKoposov et alj200^ . More recent studies have t urned their 
attention to the ultra-faint satellites: lMacci6 et alj l l2009l) have used 
three different SA galaxy formation models to study the satellite 
population of MW-like galaxies. They have used both analytic and 
numerical merging histories of MW-like haloes and shown that all 
three models reproduce the luminosity function of the MW down 
to Mv = —2, with a hint for a bending around My = —5. 

In this paper we combine high-resolution simulations of a 
MW-like halo with a S A galaxy formation model to investigate how 
various astrophysical processes affect the formation and evolution 
of satellites around the Milky Way. Ou r study extends the analysis 
presented in iDe Lucia & Helmil feOOSi) , which focused on the for- 
mation of the MW galaxy and of its stellar halo. We find that by pre- 
venting cooling in haloes with Tvir < 10* K (the atomic hydrogen 
cooling limit) and including the impact of the reionization of the 
Universe, our model is able to repr oduce the latest measu rements of 
the satellite luminosity function bv lKoposov et al.l ( l2008h . We show 
that the same model reproduces the metallicity distribution function 
(MDF) of the MW satellites, by including a route to recycle metals 
produced in newly formed stars through the hot phase. Our model 
satellites exhibit several scaling relations similar to those followed 
by the MW satellites, such as the metallicity-luminosity and the 
luminosity-size relations. The properties of the model satellites re- 
sembling the newly discovered ultra-faint SDSS satellites appear to 
be sensitive to the SN feedback recipe adopted. As we will discuss 
in the following, our model suggests that the surviving satellites are 
generally associated with haloes whose present-day peak circular 
velocity, Vinax, ^ lOkms"^, total mass exceeded a few 10® M0 
ai z ~ 10 — 20 and which acquired their maximum dark matter 
mass well above the cooling threshold, after 2; ~ 6. 

This paper is organised as follows. Section ITT] presents the 
simulations used in this study, and in Section |Z2l we summarise 
our semi-analytical galaxy formation model, emphasising the new 
features added to the model presented in lDe Lucia & Helmil ( l2008h . 
In Section[3]we present our main results and in Section|4]we discuss 
the implications of our study. We give our conclusions in Section[5] 



2 THE HYBRID MODEL OF GALAXY FORMATION 
AND EVOLUTION 

2.1 ACDM Simulations of a MW-like halo 

We h ave used a series o f high resolution simulations of a MW-like 
halo dStoehr et al]| 20021 ; [stoe hr 2006). We no te that this is the ex - 
act GAnew series used in prev ious stud ies bv lLi& Helmil ( l2008h . 
iDe Lucia & Helmil j2008l) and lLietalJ ( |2009|) . We therefore only 
summarise the basic properties of the simulations here, and re- 
fer the reader to those papers for more details. The simulations 
were carried out with GADGET-^ (Springel et al. 2001) adopting 
a ACDM cosmological model, with fio = 0.3, = 0.7, /i = 0.7, 
0-8(2 = 0) = 0.9, and Hubble constant Ho=100/i km s"^ Mpc"\ 

^ The GAnew simulations were in fact carried out using an intermediate 
version between GADGET and GADGET-2. 
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The target MW-like halo was simulated at four increasing reso- 
lution levels (the mass resolution was increased by a factor 9.33 
each time). In the highest resolution simulation (GA3new), there 
are approximately 10^ particles with mass nip — 2.063 x 10^ Mq 
h^^ within the virial radius. Each re-simulation produced 108 out- 
puts, equally spaced logarithmically in time between z = 37.6 and 
z = 2.3, and nearly linearly spaced from z — 2.3 to present. 

Merger trees for all self-bound haloes w ere constructed as de- 
scribed in detail in iDe Lucia & Helmil ilOO^ . Virialised structures 
were identified using the standard friends-of-friends (POP) algo- 
rithm and linking all particles separated by less tha n 0.2 the mean 
inter- particle separation. The algorithm SUBFIND dSpringel et al.l 
l200lh was then applied to each POP grou p to find the gravita - 
tionally self-bound substructures. Pollowing lNavarro et al.l ( fT997h . 
we define -R200 as the radius of a sphere enclosing a mass, M200, 
whose average density is 200 times the critical density of the Uni- 
verse at redshift z, i.e. : 



M200 = 



47ri?ioo 



■200Scrit{z) 



)^^200 



G 



The velocity, V200, is defined as the circular velocity of the halo at 
R200 (V200 — \/ GM2oo/R20o)- M200 is directly measured from 
the simulations and used to calculate -R200 and V200. In our mod- 
els, we approximate the virial properties of dark matter haloes, e.g. 
virial radius {Rvir), virial mass (A/vir) and virial velocity (Vvir) 
by -R200, A/200 and V200 respectively, unless otherwise explicitly 

stated. 

Pollowing lOe Lucia & Helmil ( |2008|) . we scaled down the 
original outputs by a factor of 1.42^^ for the mass and 1.42 for 
the positions and velocit ies, in order to have a MW-like halo 
with V200 ~ 15 0kms~^ jSattaglia et alj|2005l : ISmith et al.ll2007l ; 
IXue et al.ll2008h . After the scaling, the smallest resolved subhaloes 
(containing 20 particles) have dark matter mass Mum ~ 2 x 
10^ M0 in the highest resolution run (GA3new). A/dm denotes 
the total bound mass at 2: = determined by SUBFIND through- 
out this paper The present-day virial mass and the virial radius for 
the MW-like halo ai'e A/2oo ~ 10^^ M© and i?200 = 209 kpc, 
respectively. 



2.2 Semi-analytic modelling 

We use a semi-analytical galaxy formation model to study the bary- 
onic properties of a MW-like galaxy and its satellites. This model 
has been developed mainly at the Max-Planck-Institut fiir Astro- 
physik and we refer to it as the 'Munich model' later in the text. 
The essential id eas of any semi-analy tic m odel can be traced bac k 
to the works by IWhite & Reed l ll978l) and lWhite & PrenkI l ll99lh . 
and include physical processes such as the cooling of gas, star for- 
mation and feedback due to supernova explosions. Over the years, 
the Munich model has been enriched wit h new 'ingredients', e.g. 
the g rowth of supermassive black holes jKauffma nn & Haehnell 
200c ), the inclusion of dar k matter substructure s |springel e t al 



20011) , chemical enrichme nt toe Lucia et alj2004l) and AGN feed- 



back dCroton etai]|200d). The model we use in this study has 
been presented in iDe Lucia & Helmil bOOSi) and builds upon 
the model whose results have been made publicly available 
toe Lucia & Blaizo3l2007h . In order to reproduce the properties 
of the MW satellites, we have made a few modifications to the 
original model. In the fo llowing, we refer to the model used by 
iDe Lucia & Helmil ( |2008|) as the MW-model and to the fiducial one 
used in this work as the satellite-model. Here we give a brief sum- 



mary of the physical processes implemented in our models which 
are crucial to the properties of the satellites. 



2.2.1 Reionization 



Pollowing Croton et alj ( |2006|) . we make use of the results by 
iGnedinI j2000l) who simulated cosmological reionization and quan- 
tified the effe ct of ph otoionization on the baryon fraction of low- 
mass haloes. iGnedinI found that reionization reduces the baryon 
content in haloes whose mass are smaller than a particular 'filter- 
ing mass' scale, that varies with redshift. The fraction of baryons, 
/^'"'°, in a halo of mass A'/vir at redshift z, is decreased compared 
to the universal baryon fractior0 according to the ratio of the halo 
mass and the 'filtering mass', Mf'. 



/halo / T\ f 
h (z, A/v 



fb 



[H-0.26A/F(2)/A/vir 



(1) 



Por A/f(z) , we use the analytic al fitting function given in Ap- 
pendix B o f lKravtsov et"ai] | |2004 . 

In our fiducial satellite-model, we assume that reionization 
starts at redshift zo — 15 and ends at Zr = 11.5, while the original 
MW-model adopted zq = 8 and z,. — 7 (in both models reioniza- 
tion lasts for about 0.12 Gyr). Por simplicity, we refer to the reion- 
ization epoch, Zicio, as zo hereafter Our choice of fixing Zrdo = 15 
gives the right shape and the normalisation for the satellite luminos- 
ity function. We discuss the dependence of the satellites luminosity 
function on different assumptions for the reionization epoch in Sec- 
tionO 



2.2.2 Cooling 

In our model, the cooling of the shock-heated g as is treated as a 
classical cooling flow (e.g. I White & PrenkI 1 9^ ), with the cooling 
rate d epending on the temp erature and metallicity of the hot gas. 
As in iDe Lucia et al.l | |2004) . we model thes e dependences using 
the co llisional ionisation cooling curves of ISutherland & Dopital 
( Il993h . Por primordial (or low-metallicity) composition, the cool- 
ing is dominated by bremsstrahlung (free-free) emission at high 
temperatures (T > 10*^ K) and it is most efficient at T ~ 10^ K 
and ~ 1.5 x 10* K (the H and He^ peaks of the cooling func- 
tion). Line cooling from heavy elements dominates in the 10^-10^ 
K regime for non primordial compositions. Por T < IC* K, i.e. 
below the atomic hydrogen cooling limit, the dominant coolant is 
molecular hydrogen (H2). The virial temperature of a halo can be 
expressed as a function of its virial velocity as: 



Tvir(2) = 35.9 



v.iAz) 

km s^^ 



(2) 



Therefore, a halo with T^ii = IC* K corresponds to a virial ve- 
locity of Vvir = 16.7 km s"'^, which is equivalent to A/vir ^ 
3 X 10^ M0 when z = 15, and A/vir ~ 2 x 10^ Mq when z = 0. 

In the MW-model, haloes with Tvir lower than 10* K, are 
able to cool as much gas as a 10* K halo with the same metal- 
licity. In the satellite-model, we forbid cooling in small haloes 
with Tvir < 10* K, for any metallicities and at all times. This 
is a reasonable approximation since molecular hydrogen is very 
sensitive to photo-dissociation caused by UV photons from (first) 



We use the WMAP 3-year value = 0.17 jSpergel et alj2007l) . 
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stellar objects i Haiman et alj|2000l . see also iKravtsov et aLll2004l : 
IKoposov et al.ll2009t) . 



2.2.3 Star formation and supernova feedback 

T he star forma tion model used in this work i s described in detai l 
in lDe Lucia'&'Helmi (2008h . while we refer to lCroton l l20Qd) 
and De Lucia et alj ( |2004|) for details on the supernova feedback 
models. Cold gas is assumed to be distributed in an exponential 
disk and provides the raw material for star formation, which occurs 
at a rate: 



(3) 



where asp is a free parameter which controls the star formation 
efficiency, tdyn is the disk dynamical ti me and Mat is the gas mass 
abo ve a critical density t hresho ld. As in lOe Lucia & Blaizoll ^20071) 
and iDe Lucia & Helmil | |2008|) , we fix asF at 0.03, and adopt a 
Chabrier IMF. The surfa ce density threshold takes a constant value 
throughout the disk (see lOe Lucia & Helmiir2008l for details). At 
each time step. At, we calculate the amount of newly formed stars. 



AM* =jp At. 



(4) 



Massive stars explode as SNe and inject energy in the sur- 
rounding interstellar medium. In our model, we do not consider the 
delay between star formation and SN energy (and metals) injection, 
i.e. the lifetime of such stars is assumed to be zero. The energy in- 
jection by SNe per solar mass can be expressed as Vg^^ = rjsN-EsN 
where tisn — 8.0 x 10^"^ the number of SNe per unit solar 

mass expected from a Chabrier IMljf] and Esn = 1-0 x 10^^ erg 
is the average energy per SN. The energy released by SNe in the 
same time interval is: 



AEsN = ehaio ■ 0.5\/snAM. 



(5) 



where ehaio represents the efficiency with which the energy is able 
to reheat disk gas. We assume that the amount of cold gas reheated 
by SNe is proportional to the newly formed stars: 



AM,, 



sk AM. 



(6) 



If this gas were added back to the hot phase without changing its 
specific energy, the total thermal energy would change by: 



(7) 



If AEsN > Ai5hot, SN feedback is energetic enough to eject some 
of the hot gas outside the halo and we assume that: 



AM 



eject 



A^sN - Ag hot 
O.SKl 



= Ehalo 



V3 



EdiBk AM*. (8) 



The ejected material can be re-incorporated into the hot compo- 
nent associated with the central galaxy as the ha lo keeps grow- 
ing by accreting materia l from the surroundings jPe Lucia et al.l 
l2004l : ICroton et all 1 20061) . For the two parameters which regulate 
the feedback, ehaio a nd edisk, we assume t he va lues 0.35 and 3.5 re- 
spectiv ely, following ICroton et alj | |2006|) and IPe Lucia & Blaizo^ 
( |2007|) . This means that for galaxies with Vvir < 200 km s , 
AMeject > for AM. > 0. 



^ This is the fraction of stars with mass larger than ~ 8 Mq per unit of 
stellar mass formed. Note that, since we have adopted an instantaneous re- 
cycling approximation, it is more appropriate to compare the metallicities 
of our model galaxies with elements synthesised by Type II supernovae. 



An alternative supernova feedback 

Since dwarf galaxies have shallow potential wells, their properties 
are expected to be particularly sensitive to the adopted SN feedback 
model. To explore this dependency, we have also used an alternative 
SN feedback recipe in addition to the 'standard' recipe mentioned 
above. We note that this altern ative feedback recipe is equivalent to 
the ejection model described in lPe Lucia et al. Illooi), and we refer 
to this paper for more details on this particular feedback model. In 
the ejection model, the gas reheated by SNe is computed on the 
basis of energy conservation arguments and depends on the galaxy 
mass as: 



vir 



(9) 



Our choices of the feedback parameters (e — 0.05 and Vsn ~ 
634kms^^) imply that galaxies with Vvir < 87kms^^ would 
have in this model more heated mass per unit of newly formed stel- 
lar mass, compared to the standard recipe (Eq. |6ll. In this model 
the material reheated by supernova explosions in central galaxies is 
assumed to leave the halo and to be deposited in an 'ejected' com- 
ponent that can be re-incorporated into the hot gas at later times. 
The material reheated by SNe explosions in satellite galaxies is as- 
sumed to be incorporated directly in the hot component associated 
with the corresponding central galaxy. 

We will discuss later how our results for the baryonic proper- 
ties of the satellites depend on the adopted feedback scheme. 



2.2.4 Metal recycling through the hot phase 

At each time step, the masses exchanged among the four phases: 
A/hot, Mcoid, M«, Mojoct, (i.e. hot gas, cold gas, s tars and ejecta) 
are up dated as described in Section 4.7 and Fig. I in lPe Lucia et al.l 
( l2004t) . The metallicity in each phase is denoted as and is de- 
fined as the ratio between the mass in metals in each component 
(M^) and the corresponding mass (Mx) where the suffix x is hot, 
cold, star or eject. In the satellite-model, we include a route to re- 
cycle metals produced by newly formed stars throu gh the hot phase 
of a g alaxy. The equations given in Section 4.7 of |Pe Lucia et al.l 
( l2004r) modify as follows: 



Mi 



+ (1 - i?) ■ V ■ ^cold 



— ^Mcool ■ -^tiot + Mback ■ Zc]cct 



+ 5](M.e 

sat 



■ ^cold) + fzhot - Y -Tp 



M' 



Mi 



+ A/cool • Zi,ot — {1 ~ R) ■ Ip ■ Zcold 

+ (1 - i^Zhot) ■ y • ^ — A/reheat ' ^cold 



= +M 



eject 



Ztiot — A/b 



Reject . 



A constant yield Y of heavy elements is assumed to be produced 
per solar mass of gas converted into stars. The gas fraction returned 
by evolved stars is 7i = 0.43, appropriate for a Chabrier IMF. In the 
above equations, Mcooi represents the cooling rate; A/back provides 
the re-incorporation rate; i\/rcheat is the reheating rate by SNe, and 
Mojoct is the rate of mass ejected outside the halo. 

For the alternative SN feedback recipe, the reheated gas is as- 
sumed to be ejected from the cold phase directly for central galax- 
ies. In this case, the metallicity in the ejecta is updated as 
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-^ojoct — +^'^rchoat ' ^cold — A^back ' ^ejoct • 

In the MW-model, all newly produced metals returned to the 
cold phas e immediately, i.e. fzhot = 0. Hydrodynamical simu- 
lations by iMac Low & Ferrari i ll 9991) suggest that metals can be 
blown out from small galaxies with gas mass below 10^ Mq (cor- 
responding to a halo of Mvir = 3.5 x 10* M©). In our satellite 
model, we assume a simple two-state value for Fzhot to account 
for the above mass dependence: 

^ _ f 0.0 if Mvir 5: 5 X 10^° Mq 

1 0.95 otherwise. 

This means that for galaxies with a dark matter halo with virial 
mass less than 5 x 10^" Mq, 95 per cent of newly produced metals 
are deposited directly into the hot phase in this model. 

2.3 Treatment of satellite galaxies 

We follow the convention established along the development of the 
Munich model to classify galaxies according to their association 
with a distinct dark matter substructure. The galaxy associated with 
the most massive subhalo in a FOF group is referred to as 'Type 
0' or central galaxy. Other galaxies in a FOF group are usually 
referred to as satellites and are further differentiated into 'Type 1 ' 
galaxies, if their dark matter subhalo is still identified, and 'Type 
2', when their subhalo has fallen below the resolution limit of the 
simulation. 

When a galaxy becomes a satellite, the dark matter mass of the 
parent subhalo is approximated using the number of bound dark 
matter particles given by SUBFIND. The disk size is fixed at the 
value it had just before accretion. In our model, only Type central 
galaxies are allowed to accrete the material that cools from the hot 
gas associated with the parent FOF group. Satellite galaxies do not 
have hot and ejected components and the cooling is forbidden, i.e. 
A/hot = Afejoct = Afcooi = 0. When a galaxy becomes a satellite, 
its hot and ejected components are transferred to the corresponding 
components of the central galaxy. As a consequence, once matter 
leaves the cold phase of a satellite, it does not rejoin the (Type 1 
or Type 2) satellite at a later time, but it can be accreted onto the 
corresponding central galaxy. 

2.4 Model improvements and the central MW-like galaxy 

The values of the parameters that enter in the SA model were cho- 
sen so as to reproduce several observations of galaxies in the lo- 
cal Universe, in particular the local galaxy luminosi ty function and 
the mass and lu minosity of MW-like galaxies (see ^ De Lucia et al] 
l2004ICrotone t al. 2006; De Lucia & Blaizot 2007). It is therefore 
not entirely surprising that this same parameter set provides results 
that are in nice agreement with ob servational properties of our MW 
galaxy l IDe Lucia & HelmiluOOSh . In order to obtain a good agree- 
ment with the observed properties of the MW satellites, however, 
we found that we had to implement some slight modifications. 

In our satellite-mode l, we adopt the same parameter set used in 
iDe Lucia & Helmil ( |2008|) except for the reionization epoch, Zrcio, 
and fraction of metal ejected into the hot component, Fzhot ■ Table[T] 
introduce s the SA mo dels that we employ in this study. In the MW- 
model bv lDe Lucia & Helmi (^008), reionization occurs at Zrcio = 
8, haloes with Tvir < 10* K are allowed to cool at the same rate 
of a 10* K halo with the same hot gas metallicity, a 'standard' SN 
feedback recipe is adopted and all new metals are kept in the cold 
gas phase (Jzhot = 0). Our fiducial model for the satellite galaxies 



corresponds to the 'satellite-model' (the fourth row in TableO. In 
this model, reionization starts at z = 15, and 95 per cent of the new 
metals are deposited directly into the hot component in galaxies 
with A/vir < 5 X 10^" Mq. In all models with the prefix 'satellite', 
we forbid cooling in haloes with Tvir < 10* K. 

Below, we briefly discuss the dependence of model results for 
the central MW-like galaxy on our changes for Zrdo, the cooling in 
small haloes and metal recycling through hot phase. The results for 
the satellite galaxies are presented in Section[3l Table|2]summarises 
the properties of the MW-like galaxies in different SA models at 

2 = 0. 

The results correspon ding to the MW-model a re given in the 
first row (see also Fig. 2 of lOe Lucia & Helmil2008h . The only dif- 
ference between the MW-model and the satellite-model A is the 
suppression of cooling in small haloes. Comparing the results of 
these two models, we find no significant changes in the properties 
of the present-day MW-like galaxy. In satellite-model B, we also 
change Zrcio to 15 and keep Fzhot = 0. The only significant effect 
of an early reionization is to bring down the black hole mass by 
~ 15 per cent compared to the value obtained in the MW-model. 
The early reionization also results in a slight increase of the stel- 
lar mass but this is still well within the observational uncertainties 
(M. ~ 5 - 8 X 10^° Mq). 

Our fiducial satellite-model gives a total stellar mass similar 
to that obtained from the MW-model and in agreement with cur- 
rent observational constraints. The results of these two models are 
also very close in terms of the mass of the bulge and the cold gas 
content. The black hole mass Mbh ~ 6.9 x 10^ Mq from the 
satellite-model is in marginal agreement with the latest measure- 
ment of the MW black hole mass A/bh = (4.5 ± 0.4) x lO'' Mq 
iGhez et al.|2008| ). The ejection of metals into the hot component in 
small galaxies (cf. satellite-model B) only makes the bulge slightly 
more metal-poor (~ 0.06 dex). Results listed in Table[2]show that 
the modifications discussed above influence only the properties of 
dwarf galaxies, whil e preserving the p r operti es of the MW-like 
galaxies discussed in iDe Lucia & Helmil ( |2008|) . 

The 5th row in Tables [T] and |2| corresponds to a model which 
incorporates the 'alternative' (or ejection) SN feedback scheme de- 
scribed in Section l2.2.3l When compared to the 'standard' scheme, 
the dependency of the amount of reheated gas on 1/Vjir in this 
scheme results in a more efficient feedback for small galaxies and 
in a less efficient ejection for more massive systems like the MW 
galaxy. These galaxies tend to have larger stellar masses, a more 
massive bulge and tend to be more metal-rich in this scheme. As we 
will show in Section |3] however, this alternative feedback scheme 
is able to better match the properties of the MW satellites. 

We therefore propose a combination of these two feedback 
recipes to account for the properties of galaxies on large and small 
mass scales. This model corresponds to the 'satellite-model com- 
bined', and in this scheme, we calculate the amount of gas reheated 
by SNe depending on the local potential well (mass of the associ- 
ated subhalo): 

I EdiskAA^, Otherwise. 

The reheated gas is treated as in the ejection scenario, i.e. it is added 
to the ejected component of a central galaxy and lost into the hot 
component for a satellite. The 'satellite-model combined' entry in 
Table |2| lists the properties of the MW-like galaxy in this scheme. 
As expected, the stellar mass, the total and bulge metallicity, as well 
as the total luminosity are now very similar to what we get with the 
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Table 1. Nomenclature and features of the semi-analytic models used in this study. 



Model Name 


20, 


cooling for Vvir < 16.7 km s ^ 


SN feedback 




(1) 


(2) 


(3) 


(4) 


(5) 


MW-model 


(8. 7) 


Yes 


Standard 


0.0 


satellite-model A 


(8,7) 


No 


Standard 


0.0 


satellite-model B 


(15, 11.5) 


No 


Standard 


0.0 


satellite-model 


(15, 11.5) 


No 


Standard 


(0.95, 0.0) 


satellite-model ejection 


(15, 11.5) 


No 


ejection 


(0.95, 0.0) 


satellite-model combined 


(15, 11.5) 


No 


combined 


(0.95, 0.0) 



Different columns list: (1) the model name; (2) the reionization epoch; (3) the adopted cooling recipe in haloes with Tvir < 10^ K; (4) the adopted feedback 

recipe; (5) the fraction of metals injected directly into the hot component. 



Table 2. Properties of the Milky Way-like galaxy in the semi-analytic models employed in this study, for the GA3new simulation. 



Model Name 
(1) 


M, 
[IQio M©] 
(2) 


A'/bulge 

[IQIOMq] 
(3) 


-^^coldgas 

[IQIOMq] 
(4) 


Mbh 
[10^ Mq] 
(5) 


[dex] 
(6) 


[dex] 
(V) 


Mb 
[mag] 
(8) 


MW-model 


5.73 


0.64 


1.06 


8.2 


-0.05 


-0.28 


-19.53 


satellite-model A 
satellite-model B 


5.75 
5.85 


0.64 
0.62 


1.07 
1.11 


8.0 
7.1 


-0.05 
-0.06 


-0.28 
-0.30 


-19.50 
-19.52 


satellite-model 


5.88 


0.63 


1.11 


6.9 


-0.06 


-0.35 


-19.52 


satellite-model ejection 


8.26 


2.58 


1.09 


12.8 


0.12 


0.00 


-19.08 


satellite-model combined 


5.02 


0.84 


0.95 


12.0 


-0.08 


-0.25 


-19.53 



Different columns list: (1) the model name; (2) stellar mass; (3) mass of the bulge; (4) cold gas content; (5) mass of the black hole; (6) logarithmic value of 
the total stellar metallicity; (7) logarithmic value of the bulge metallicity; (8) B-band absolute magnitude connected for internal dust attenuation. 



standard feedback scheme, albeit the bulge (and the black hole) are 
now more massive. The last two models listed in Table [T]populate 
the same set of subhaloes with stars, with almost identical proper- 
ties. For simplicity, we will not discuss the results of the combined 
scheme for satellite galaxies. 



more than 97 per cent of their dark matter by present. Therefore 
it is quite unlikely that a bound stellar core would survive such 
a severe tidal stripping. We will discuss more about Type 2 model 
galaxies in Section |4~2l Throughout this paper, we will refer to sub- 
haloes that host stars as 'luminous satellites' or simply 'satellites', 
and will refer to those that do not host any star as dark satellites or 
subhaloes. 



3 THE IVIILKY WAY SATELLITES 

In this Section, we define as SA model satellites of the MW-like 
galaxy those that satisfy the following conditions at 2 = 0: (i) 
a satellite has to belong to the same FOF group where the IVIW- 
like galaxy is; (ii) the distance to the MW-like galaxy must be 
< 280 kpc; (iii) the galaxy is associated with a dark matter sub- 
halo (i.e. it has to be Type 1 galaxies). The distance cut corre- 
sponds to the current observational limits, but we include in our 
comparison the very distant and recently discovered satellite Leo T 
(at ~ 420 kpc from the Milky Way). 

We do not consider here satellites that had their dark matter 
haloes tidally stripped below the resolution limit of the simulation 
(Type 2 galaxies). This selection is motivated by the fact that while 
resolved, our Type 1 galaxies are dark matter dominated at all radii 
(see Section [33] and in agree ment with observations of satellites 
around the MW and M31, e.g. Isimon & Gehall2007l ; IStrigari et all 
l2007ll2008l :l Walker et al.l2009h . Since Type 2 galaxies in our fidu- 
cial model, reach a maximum virial mass during their evolution 
larger than 6.8 x 10^ M©, this means that these galaxies have lost 



3.1 The satellite luminosity function 

Our fiducial satellite-model gives 51 luminous satellites within 
280 kpc for GA3new. This is in good agreement with the estimated 
'al l sky' number of sate llites (--^ 45) brighter than My ~ —5.0 
bv lKoDosov et alj ( l2008h . If we remove the distance constraint, the 
number of satellites is only increased by one, and the number of 
subhalos varies from 1865 to 1869. 

The mass functions for the fiducial satellite-model and of the 
surviving subhaloes within 280 kpc in GA3new are shown in Fig.[T] 
The mass plotted here is the dark matter mass at 2 = determined 
by SUBFIND. As indicated by the dashed histogram, all subhaloes 
with present-day Mdm > 10^ Mq resolved in GA3new host lumi- 
nous satellites. The mass function of these subhaloes deviates from 
the power-law shape mass function of the full subhalo population 
and is fairly flat below i\/DM = 10^ M©, down to the resolution 
limit (Mdm ~ 10^'^ M©). For comparison, the dotted histogram 
in Fig.[T]shows the mass function of surviving subhaloes, within the 
same distance range from the central galaxy, from the lower reso- 
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Figure 1. The solid histogram shows the present-day mass function of 
model satellites using the satelUte-model for the highest resolution sim- 
ulation. Error bars denote the l-cr Poisson uncertainties. The dashed his- 
togram shows the subhalo mass function for the highest resolution simula- 
tion, which steeply rises up to the resolution limit. The dotted histogram is 
the subhalo mass function for a lower resolution simulation (i.e. GA2new). 



Table 3. Number of satellites around the model MW galaxy for the different 
SA models used in this study for the GA3new simulation. 



Model Name 




20, Zr 




(1) 


(2) 


(3) 


(4) 


MW-model 


286 


(8,7) 


0.0 


satellite-model A 


88 


(8, 7) 


0.0 


satellite-model B 


51 


(15, 11.5) 


0.0 


satellite-model 


51 


(15, 11.5) 


(0.95, 0.0) 


satellite-model ejection 


51 


(15, 11.5) 


(0.95, 0.0) 


satellite-model combined 


51 


(15, 11.5) 


(0.95, 0.0) 



Different columns list: (1) the model name; (2) the number of luminous 
satellites; (3) the adopted reionization epoch; (4) the fraction of metals 
ejected directly into the hot component. 

lution simulation GA2new. The smallest subhalo which could be 
resolved in GA2new has A/dm ~ 2 x 10^ Mq . The subhalo mass 
functions from the two simulations agree well down to 10* M©. 
At lower masses, numerical effects start to become important for 
GA2new. However, the number of satellites with M^m < 10^ Mq 
is still lower than the number of subhaloes resolved in GA2new, 
which suggests that numerical resolution should not be an issue 
and that the observed decline of luminous satellites is a result of 
how we model the baryonic physics, e.g. SN feedback, as we will 
see later. 

Fig. [2] compares the mass function for the SA models ex- 
plored in this study. In Table [3] we list the number of luminous 
satellites for these models. Note that all the models which adopt 
Zroio = 15 and forbid cooling in haloes with VvIy < 16.7kms~^, 
i.e. satellite-model B, satellite-model, satellite-model ejection and 
satellite-model combined, populate galaxies in the same set of 51 
subhaloes. The mass functions of these four models are therefore 
identical. 

Fig. [3] shows the luminosity functions of the first three SA 
models listed in the tables. The (red) dashed line shows the 'all sky 



Figure 2. Present-day mass functions (solid histograms) for dark matter 
subhaloes associated with satellites predicted by the semi-analytic models 
listed in Table [T] The dashed histogram again indicates the subhalo mass 
function for the GA3new simulation. 

SDSS' powe r-law luminosity func tion of satellites within 280 kpc, 
estimated bv lKoposov et al.l ( l2008h . For satellite-model B, we show 
the l-(j Poisson noise. The filled circles show the luminosity func- 
tion of the 22 known satellites of the MW, in cluding the latest ultra- 
faint satellite Leo V teelokurov etai]|2008l) . As a reference, in the 
MW-model the number of surviving satellites is 286. The drastic 
difference between the MW- and the satellite-model B is mostly at 
the faint end of the luminosity function, and it is due to the com- 
bined effect of an early reionization and no cooling in small haloes. 
In a model with no cooling in small haloes and a later reioniza- 
tion, the number of satellites is 88 (compare the MW-model and the 
satellite-model A). The number of satellites is further reduced to 5 1 
when assuming an earlier reionization epoch (zrcio = 15), with a 
reduction of galaxies in the luminosity range Mv G [—7,-10] 
compared to a model with a later reionization. We have also exper- 
imented values for 2i-cio at redshift 10 and 12 and found a number 
of surviving satellites of A^sat = 79 and 58, respectively. To sum- 
marise, these choices of z-^do ~ [8, 10, 12, 15] all give a number 
of satellites down to My = —5 which is consistent with the es- 
timation by iTollerud et al.l l l2008h (see their Fig. 6). However, af- 
ter examining the shape of the luminosity functions with differ- 
ent Zrcio, we decided to use Zi-cio = 15 in our fiducial model, as 
this choice gives both the right normalisation and a shape which is 
in better agreement with the observational measurements (see next 
paragraphs). We recall that our results are based on only one halo, 
with mass comparable to that estimated for our Milky Way. A cer- 
tain scatter in model predictions, due e.g. to the assembly history 
of the parent halo, is expected and this might be significantly larger 
than the Poisson noise plotted in Fig. [3] 

The left panel in Fig. |4] shows the luminosity function of the 
satellites in our fiducial satellite-model (solid hist ograms), com- 
pared to the 'all sky SDSS' luminosity function bv lKoposov et al.l 
1 I2OO8I) and the observational data for the 22 known MW satellites. 
The model luminosity function covers a similar luminosity range 
as the 22 MW satellites though the model predicts a lower num- 
ber of faint (Mv > — 5) satellites with respect to the expected 'all 
sky' luminosity function. In fact, the satellite-model does not pre- 
dict any satellite fainter than i\fy = —4 mag. On the other hand, 
it shows an excess with respect to the data of 10 — 15 satellites 
around Mv ~ —10. These satellites in excess are all within a dis- 
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Figure 4. Luminosity function for the observed MW satellites (filled circles and red dashed fine) and for the model satellites (black histograms). The left panel 
is for our fiducial model and the right panel is for the model with a more efficient SN feedback for small galaxies (see Section [2.2.3t . Data ai'e the same as 
those in Fig.|3] 
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Figure 3. Luminosity function of the observed MW satellites and for the 
model satellites around a MW-like galaxy, for the first three models listed 
in Table[3] The integrated y-band luminosities of MW sate llites (dots) ar e 
taken from various sources. The classi cal dSphs are from [Mated jl998l) ; 
most of the ultra-faint dwarfs are fro mlMartin et alj <2008l l. except Leo T 
<Rvan- Weber et al.l2008ll and Leo V teelokurov et alj2008h . 



tance of 200 kpc from the host galaxy and have half-light radii be- 
tween 100 and 500 pc. They should therefore have appeared in the 
'all-sky' SDSS luminosity function if they existed (provided the 
distribution for these bright satellites is isotropic, as assumed by 
IKoposov et al.ll2008l) . It also seems that the model under-predicts 
the number of very bright satellites. This is, however, the regime 
where both the data and the simulations suffer from small number 
statistics. 

A comparison to Fig. [5] shows that the luminosity function 
given by satellite-model B and by satellite-model are very simi- 
lar, implying that satellite luminosities do not depend strongly on 
the fraction of newly produced metals put into the hot component 
of a galaxy. However, the metallicity distribution functions differ 
significantly (see Section [T2t . 

The right panel of Fig. |4] shows the luminosity function of 



surviving satellites resulting from the alternative feedback scheme 
(satellite model ejection). This model luminosity function agrees 
very well with the observations. When compared with that from 
the standard feedback scheme, the alternative luminosity function 
extends to fainter luminosities, reaching Mv ~ —3 and does not 
show any excess at My ~ —10. Note that in the satellite-model 
ejection, the same 5 1 subhaloes are populated with luminous galax- 
ies. This suggests that the SN feedback alone is unlikely to solve the 
'missing satellites problem' l lSomervillell2003) . and that the pres- 
ence/absence of a luminous galaxy within a dark matter substruc- 
ture is due to the particular assembly and dynamical history of the 
halo, and to the reionization history of the Universe. 

These different results are entirely due to the SN scheme 
adopted. In the ejection scheme, the amount of gas reheated by 
SNe scales as l/V^^j^, and for a galaxy with Vvir < 87kms~^, 
more gas is heated by per unit of newly formed stars with respect 
to the standard scheme, where the amount of gas heated by SNe is 
only proportional to the amount of newly formed stars. In this latter 
case, for galaxies with low star formation rates (as is the case for 
galaxies that live in small subhaloes), only very little gas is heated. 
As long as the density threshold for star formation is met, these 
galaxies keep forming stars. This is most likely the reason why our 
fiducial model exhibits an excess of satellites around My ~ — 10 
and a lack of ultra faint objects in the luminosity function. In the 
alternative feedback model, when star formation occurs in a cen- 
tral galaxy, some cold gas is ejected outside the halo, and has to 
wait several dynamical time-scales to be reincorporated into the 
hot halo, delaying any subsequent star formation. Furthermore, the 
impact on the star formation of a satellite galaxy associated with 
a subhalo with Vvir < 87kms^^ is more drastic in this model, 
because the ejected gas cannot be re-incorporated onto the satel- 
lite any longer Small satellites therefore consume their cold gas 
reservoir more efficiently in this ejection scenario, which causes 
the surface densities of the cold gas to fall below the criterion for 
forming stars. 

Given the reasonable match of the satellite-model and of 
satellite-model ejection to the observed luminosity functions, in the 
rest of the paper we will concentrate on the properties of the satel- 
lite populations in these two SA models. 
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Figure 5. Histogram of the mean iron abundance [Fc/H] determined for red giant branch stars in the MW satelhtes. For the model satellites we plot 
log(Z»/ Zq). The left panel compares the MW satellites with model satellites from our fiducial model and the right panel from the model with a more 
efficie nt SN feedback for dwarf galaxies (see Section |2. 2. 3[. Data for the MW satellites are taken from various sources: LMC and SMC from Westerlund 
il997l) : Sgr from^Cole (20 01J; Ursa M inor and Dra co fromlHarbeck et al ll l200ll) : Sextans, Sculptors, Carina and Fornax from the DART survey iHelmi et alj 
200d): Leo 11 from iKoch et al.l j2007h : Leo I from lKoch et alj j 2007 ) . For the newly disco vered SDSS ultra-f aint dwaifs, we take the measureme nts from 
Kirbv et alj j2008h . except for Bootes I for which we use lMufioz et alj J200^ . Bootes II from lKoch et alj j2009l) and Leo V from lWalker et alj j2009l) . 



3.2 The metallicity distribution 

The left panel of Fig. |5] compares the metalhcity distribution of 
model satellites in our fiducial satellite-model with the observed 
distribution. The metallicity in our model is mass-weighted and is 
defined as the ratio between the mass of metals in stars and total 
stellar mass: 

Since our model does not distinguish the long-lived main iron con- 
tributors (SN la ) from the short-lived a-elements enrichers (SN II ), 
a direct comparison of Z to [Fe/H] is questionable. Nevertheless, 
here we assume that the logarithmic value of the mass-weighted 
metallicity normalised to the solar value (Zq = 0.02) can be com- 
pared qualitatively with the [Fe/H] derived from spectra of Red 
Giant Branch stars in the MW satellites. 

Among the 5 1 surviving satellites in the satellite-model, four 
of them are free of metals since they have only made stars once 
from pristine gas. These four metal-free satellites all have present- 
day A/dm < 10'-* Mq and M, of lO" - lO'"^ Mq. We do not in- 
clude these 'metal-free' satellites in the left panel of Fig. [5] and 
(metallicity) related discussions. Fig. |5] shows the histograms of 
the mean [Fe/H] of resolved stars in each MW satellite. The dis- 
tribution of the 1 1 classical MW satellites is shown by the dotted- 
dashed histogram, while the corresponding distribution including 
also 1 1 ultra-faint satellites is given by the dashed histogram. The 
metallicity distributions of model and MW satellites cover similar 
ranges. However, the peak of the metallicity distribution for the 22 
MW satellites is shifted to lower values with respect to the corre- 
sponding distribution from the satellite-model. The excess of model 
satellites in the range of —2 < log(Z,/ Zq) < —0.5 conesponds 
to the bump at My ~ — 10 seen in the luminosity function. 

If we had set fzhot = (as in satellite-model B), the pre- 
dicted metallicity distribution would be shifted towards even higher 
metallicity. In this case, 36 satellites are more metal-rich than 
log(.Z*/Z0) = -1, and only 11 of them have log(Z*/ZQ) < 
— 1, inconsistent with the observations. This is why our fiducial 
model is preferred to satellite-model B. 



In the right panel of Fig. |5] we plot the metallicity function 
of the 38 non metal-fre^ satellites in the satellite-model ejection. 
The peak of this distribution is also shifted to slightly higher values 
than observed, but it has a more even distribution compared to the 
standard recipe. As discussed in Section [XT] this is due to the fact 
that satellites with Vvir < 87kms^^ loose more of their cold gas 
reservoir (and a larger fraction of their metals, see Section [2.2.4t in 
this model with respect to the standard feedback scheme. 

3.3 Star formation histories 

In Fig.|6]we present the evolution of the stellar mass and the SUB- 
FIND dark matter mass predicted in our fiducial satellite-model. 
We remind that this model does not give any satellite fainter than 
AIv = —4, and gives only 2 satellites fainter than Mv = —5. We 
therefore restrict our comparisons to the classical MW satellites, 
and we sort them by their present-day luminosity into three bins: 
luminous (—16 < My < —13, similar to Sagittarius and Fornax), 
intermediate (—12 < My < — 10 like Leo I and Sculptor), and 
faint (-10 < Mv < -8 i.e. Leo II, Sextans, Carina, UMi and 
Draco). Satellites belonging to these three bins are shown from top 
to bottom in Fig.|6] We do not include here the two model satellites 
(one in the luminous and one in the intermediate luminosity bin) 
which have more stellar mass than dark matter mass at 2 = 0, due 
to significant tidal interactions with the main halo. We recall that in 
our modelling, we do not account for tidal stripping of stars and for 
the loss of the cold gas due to e.g. ram-pressure. Therefore, we also 
refrain from considering systems which form > 50 per cent of their 
stars after becoming satellites of the MW-like halo. This choice is 
motivated b y the prevalence of old stars (i.e. > 10 Gyr) seen in the 
MW dSphs jPolphin et al.l2005l : IOrban et al.l2008l) . When exclud- 
ing these galaxies, the number of satellites is reduced from 16 to 
13 in the intermediate bin; 14 to 7 in the low luminosity bin. We 
will later see that faint satellites are mostly accreted before z = 1 
thus their gas should have been reduced due to the interactions with 

^ i.e. 13 satellites are 'metal-free' in this case. 
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Figure 6. Evolution of the stellar mass (left) and dark matter mass (right) for satellites in the fiducial satellite-model. Satellites are sorted by Mv{z = 0) into 
three groups: —16 < My < —13 (top panel); —12 < My < —10 (middle panel); —10 < My < —8 (bottom panel). Stellar masses are normalised by 
the present-day values, M* (2 = 0). Different colours coiTespond to different satelHtes, and the arrows indicate the accretion time, defined as when a satellite 
was identified as a central galaxy for the last time. The same colour is used to plot the stellar and the dark matter mass for a given satellite in the panel aside. 
The vertical dashed lines mark the end of the reionization, Zr = 11.5, in the satellite-model. 



the central galaxy and the star formation should have (on average) 
ceased shortly after being accreted. 

It is encouraging that the model satellites show in Fig. |6] all 
contain stars older than 10 Gyr regardless of their luminosities, 
in good agreement with observation^. The most luminous satel- 



lites build up their stellar content over a longer period of time com- 
pared to the faintest ones as observed in the Local Group satellites 
(Dolphin et al. 20o3)- The five systems that have been accreted ear- 
liest (> 9 Gyr), are associated with fainter satellites (Aly > —12), 
and are dominated by old stars, which means that these galaxies 



^ In the full sample of model satellites, 43 out of 5 1 made their first gen- populations are excluded in the analysis here with the criterion that half of 

eration of stars 10 Gyr ago. Those that are not dominated by old stellar the stars were in place before the accretion. 
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Figure 7. Normalised cumulative radial distributions of the MW satellites 
(dashed line) compared with that for model satellites (solid line) and for the 
dark matter subhaloes (dotted line) at 2 = 0. The distance from the Sun to 
the Galactic centre is taken as Rq = 8.0 kpc. 



Stopped forming stars soon after the accretion. The most luminous 
model satelhtes are associated with massive dark matter subhaloes 
(Mdm > 10'' M0) at present time. The Figure also shows that all 
bright satellites consist of less than 1 per cent of stars made by the 
end of the reionization {z = 11.5 equivalent to a lookback time 
of 13 Gyr), while a few of the faint ones contain > 50 per cent 
of such very old stars. None of the most luminous satellites was 
accreted before z ~ 1, and they have the most extended star for- 
mation histories and reside in the most massive haloes. This is a 
natural consequence of the fact that subhaloes that were accreted 
early, have typically experienced a more significant mass l oss than 
those accreted later jPe Lucia et alJl2004[Gao et al.ll2004l) . As ex- 
plained earlier, no cooling is allowed on satellite galaxies in our 
model. After accretion, these galaxies keep forming stars, with the 
size of their gaseous disk fixed at the value when they were ac- 
creted. The surface density of the gas soon drops below the density 
threshold for star formation, therefore rapidly quenching any fur- 
ther star formation in these galaxies. 



3.4 Other properties of the satellites 

3.4.1 Radial distribution 

Fig. |7] shows the cumulative radial distribution of the MW satel- 
lites as a dashed line, and of the model satellites as a solid line. 
Here we only compare our model results with satellites observed 
within a galactocentric distance of 280 kpc (i.e. excluding Leo T). 
The median distances of the two distributions agree well. For com- 
parison, we also show the distribution of all surviving dark matter 
subhaloes (with or without stars) as a dotted line. It is clear that 
the radial distribution of subhaloes is much less concentrated com- 
pared to that of the MW and the model satellites. The reason why 
our model satellites show a cumulative radial distribution similar 
to that of the MW satellites can be understood from Fig. [S] This 
Figure shows the present-day Mdm as a function of the accretion 
time. The small grey circles are for the 'dark' satellites, and the 
larger black symbols are for those subhaloes which host luminous 
satellites. The bulk of dark satellites (i.e. subhaloes without stars) in 
the mass range of 10^ M© < Mdm < 10^ M© have been accreted 
in the last two Gyr, and are preferentially found in the outskirts of 



Figure 8. Present-day bound dark matter mass as a function of the accretion 
time. Black symbols represent the luminous model satellites, while grey 
symbols mark dark matter subhaloes which failed to form stars. 



the MW halo as, for this mass scale, dynamical friction is not im- 
portant. Since these small dark satellites dominate by number, this 
leads to a much more even distribution for subhaloes than for the 
luminous satellites. Half of the satellites were accreted more than 
7 Gyr ago, and the most massive ones (e.g. A/dm > 10® M©) at 
z = Q were all accreted after z ^ 1 (see also Fig.|6j. There is also a 
clear bias for the least massive subhaloes hosting stars to have been 
accreted earlier (as mentioned in Section[33t. 



3.4.2 Luminosity-size relation 

Here we compare the half-light radii, ri/2, of model satellites with 
the observed distribution as a function of the total V^-band absolute 
magnitude. The left panel in Fig.|9]is for our fiducial model, while 
the right panel is for the alternative feedback scheme. The size of a 
satellite in our models is obtained assuming the stars are distributed 
in an exponential disk and is roughly proportional to the virial ra- 
dius of the associated dark matter subhalo, i.e. ri/2 ~ 1.2Ai?2ooEl 
where A is the spin parameter of the associate DM halo. In Fig. |9] 
we show the classical dSphs and the ultra-faint satellites with as- 
terisks and diamonds respectively. The black circles show model 
satellites. The half -light radii of the u ltra-faint satellites are from 
the recent study bv lMartin et al.l ( l200^ based on SPSS d ata, while 
those of the classical ones are from van den Berghl ( l2000l) . 

Model satellites with either of the feedback recipes have sizes 
between ~ 30 to 2000 pc, in quite good agreement with the MW 
satellites, especially the classical ones. For the luminous satellites, 
both models extend to larger sizes. The lack of faint satellites be- 
low A4v = —5 in the fiducial model makes the comparison with 
the ultra-faint satellites inconclusive. The compact satellites with 
ri/2 < 100 pc are a few magnitudes brighter than the ultra-faint 
MW satellites of similar size in this model. 

If we compare the sizes of galaxies associated with the same 
DM subhalo in the two feedback schemes, those given by the alter- 
native feedback recipe (right panel) are a bit larger than in the stan- 
dard feedback model. This is expected as in the alternative feedback 
scheme, more gas is ejected by SNe in small (central) galaxies and 



^ ^1/2 ~ f-68 rj3 for an exponential disk, and is the scale length of 

the disk which is computed as ~ -7=i?200 to the first order. 
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Figure 9. Luminosity as a function of the half-light ra dius for the MW and model satellites. Observational measurements are taken from various sources: For 
the classical dSphs, data are from the com pilation of lvan den Berghl (20003. Data for the ultra-faint dwaifs are from lMartin et al, t200S) . except for Leo V 
which is taken from lBelokurov et alj j2008l) . See the caption of Fig.|4]for the sources of the total V-band luminosity. 



this is reincorporated later, after the halo has grown in mass and 
size. Since the alternative feedback recipe pushes the luminosity 
function towards the fainter end, the comparison with the ultra- 
faint satellites is in better agreement than in the standard recipe. 
In addition, there are some faint satellites that are larger than the 
observed ultra-faint systems with central surface brightness fainter 



than the current observational limit /i„ 
dotted lines in the Figure). 



30 mag arcsec (see 



3.4.3 Metallicity -luminosity relation 

The classical MW satellites are known to follow a metallicity- 
luminosit y relation wit h metallicity increasing with increasing lu- 
minosity l lMateo|[l998l) . Fig. 1101 shows the distribution of model 
satellites (circles) in the metallicity-luminosity plane, together with 
the observational measurements for classical (asterisks) and new 
(diamonds) MW satellites. Black circles are for the fiducial model 
with the standard feedback scheme, and grey circles are used for 
the ejection model with the alternative feedback scheme. Again, we 
compare Z,/ Z0 for the model satellites with the averaged [Fe/H] 
for stars in the MW satellites so care is needed when interpret- 
ing these results. The error bars on the [Fo/H] values indicate the 
spread of the distributions in each MW satellite rather than mea- 
surement uncertainties. The model satellites in the fiducial model 
follow a similar trend as the classical MW satellites, while it is not 
possible to compare results from this model to the bulk of ultra- 
faint satellites since the model does not predict galaxies with lumi- 
nosity less than Lv ~ lO'*'^ Lq (see also Fig. |4]l. The excess of 
model satellites seen in the luminosity function (Fig. |4] left panel) 
around My ~ —10 corresponds to objects with log(Z»/ Z©) in 
[—2.0, —1.0] dex, and also produces a 'bump' in the metallicity 
distribution. 

With the alternative SN feedback recipe, model satellites also 
follow the relation of the classical satellites with some hints for 
metallicities larger than observed for Lv > 10*" L©. At the faint 
luminosity end, there is better agreement with the ultra-faint satel- 
lites, since the alternative feedback recipe predicts fainter and 
metal-poor satellites below Lv = 10^ Lq . From our model re- 
sults, it is not clear whether the luminosity-metallicity relation has 
a lower limit since we are affected by the resolution at the very 
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Figure 10. Metallicity-luminosity relation for model and MW satellites. 
Black circles show model satellites in the fiducial model, and grey cir- 
cles are used for the ejection model. For the model satellites, the logarith- 
mic values of the total metallicity in stars normalised to the solar value 
( Zq = 0.02) are plotted. For the MW satellites, shown in (blue) asterisks 
and diamonds, the mean [Fe/H] is used, and eiTor bars denote the disper- 
sions of [Fe/H] in each galaxy. See the captions of Fig.|3]and Fig.|5]for the 
data sources. 



low mass end. It is also not clear whether the predicted relation 
deviates from the ultra-faint satellites below Lv = 10^ L©. It 
is interesting to note that if we follow the ridge of the observed 
luminosity-metallicity relation from the high and metal-rich end 
down to Lv ~ 10"" Lq, the corresponding [Fe/H] is ~ —2.5 
which happens to be the lower limit for the empirical relation to 
convert the equivalent width of Ca II tripl et to [Fe/H] for RGB 
stars in dSph galaxies jBattaglia et al.ll2008i) . 



3.4.4 Cold gas content 

The majo rity of the MW dSph sate llites (including the ultra-faint 
ones, see iGrcevich & Putmarj|2009h are gas deficient. In Fig. [TT] 
we show the cold gas content as a function of the luminosity for 
our model satellites and compare model results with observational 
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Figure 11. Cold gas content as a function of \/-band absolute magnitude. 
Black circles are used for the standard feedback sch eme, and grey circles 
for the alternative scheme. HI masses are take n from lMateol i 19981) for the 
classical dSphs. and from iBriins et alj j2005h for the Magellanic Clouds. 
Apart from the gas rich Magellanic Clouds and Sculptor, most of classical 
dSphs are only constrained with upper limits (indicated by the downward 
pointing arrows). The distant Leo T is the only one with a clear detection of 
Mhi = 2.6 X lO'^Mp) among the newly discovered ultra-faint satellites 
iRvan- Weber et al ] |2008l) . The other ultra-faint satellites have no detected 
HI or an upper limit of Mh i < 10^ iOrcevich & Putmanlliooil) and 
are not plotted here. 



Figure 12. Vmax as a function of the y-band integrated magnitude for 
the non tidally disturbed model satellites. Black symbols show predictions 
for the fiducial model, and grey symbols are used f or the ejection mode l. 
Data for eight classical MW dSphs are taken from IWalker et al ] l l2009h . 
The lower limits on Vinax are indicated with the upward arrows except for 
Fornax dSph whose Vmax is better constrained. 



away from the giant disks, where environment is expected to play 
a less important role). 



measurements. The Figure shows that model satellites (from either 
feedback model) are much more gas-rich than the observed coun- 
terparts, by factors of a few hundred, except at the brightest end. 
The satellites from the alternative feedback recipe populate a simi- 
lar region as those from the standard recipe but there are some with 
significantly lower (only ~ 10^~* Mq) cold gas. We have verified 
that altering the reionization epoch or star formation efficiency does 
not have a significant impact on the present-day cold gas mass of 
satellite galaxies. This means that most of the cold gas in satellites 
was already in place when they were central galaxies, and that this 
gas has not been affected much by star formation and/or feedback. 
We recall that in our models, star formation only occurs when the 
mean surface density of the cold gas is above a certain threshold. 
Strong supernova feedback can in principle reduce the gas content 
in the cold phase, but most of the satellites form stars at very low 
rates, so SN feedback will be less important. In addition, increasing 
the supernova efficiency would reduce the luminosity (and metal- 
licity) of the model galaxies, ruining the agreement with observa- 
tional data shown earlier in this study. 

It is known that Local Group galaxies show a morphological 
segregation in the sense that gas-deficient dSphs are closer to the 
giant spirals (the MW and M31), while the gas-rich dirrs are more 
evenly distrib uted jMateo 199^. Recent studies, based on deeper 
HI surveys bv lGrcevich & PutmaiJ ( |2009|) also show that there is a 
clear correlation between galactocentric distance and HI content for 
dwarf galaxies in the Local Group. These are interpreted as indica- 
tions that the Local Group environment influences the gas content 
of s atellites through tidal stripping and ram -pressure stripping (see 
e.g. lGrebel et al.ll2003l : lMaver et alj|2006h . The excessive cold gas 
associated with our model satellites might be due to our neglecting 
of gas removal mechanisms, or to a simplified calculation of the 
surface density threshold (particularly for those satellites which are 



3.5 Dark matter halo mass and the dynamical properties 

As in lLi et al ] (l2009h . we trace the evolution of Vmax (i.e. the peak 
circular velocity of the associated dark matter subhalo). A/dm and 
the directly measured Mo.e in the last 2 Gyrs for each satellite. 
16 out of 51 model satellites are classified as tidally disturbed and 
these include those two model satellites which have experienced 
severe tidal stripping, i.e. M* > Mdm at z = using our fidu- 
cial model. Fig. [12] shows Vmax as a function of My for those 
35 model satellites that do not show signs of tidal disturbance. 
Black circles are for the standard feedback scheme, while grey ones 
are for the alternative feedback scheme. The latest constraints on 
Vmax for 8 classical MW dSphs based on Jeans/MCMC analysis by 
IWalker et al.l ( |2009|) are shown as asterisks. Under the assumption 
of constant velocity anisotropy, their data could constrain Vnax for 
Fornax, while only lower limits were derived for the other dSphs 
(upwards arrows in the Figure). Assuming that classical MW satel- 
lites are not tidally disturbed, our models predict that they are asso- 
ciated with dark matter haloes with Vmax ^ 10 km s^^, regardless 
the choice of SN feedback recipe. Thi s is fairly consistent with the 
constraints bv lPenarrubia et alj ( l2008l) and lWalker et alj ( |2009|) . 

Fig.ll3lshows the mass-to-light ratios calculated using masses 
within 0.6 kpc, Mo.a/ L, as a function of absolute magnitude. The 
luminosities shown are from our fiducial model. The data points for 
the eight classical M W dSph (i.e. excluding Sgr dSph) measured by 
IStrigari et al] ( l2007h are overplotted as (blue) asterisks. The lower 
and middle dashed lines correspond to constant mass values of 
Afo.6 = 6 X 10*^ Mq and Mo.e = 7 x 10^ Mq and indicate the 
upper and lower limits for the observed Mo.e - For each satellite in 
our model, two methods can be applied to measure the Afo.e: 1) to 
directly sum up bound dark matter particles within 0.6 kpc from the 
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Figure 13. Mass-to-light ratio of 8 Milky Way classical dSphs and model 
satellites as a function of luminosity. The mass corresponds to the dark 
matter mass within 0.6 kpc. For the MW dSphs, data are taken from 
IStrigari et all J2007h and plotted as asterisks. The lower and middle dashed 
lines correspond to constant values of A/o.e = 6 X 10® Mq and 7 X 
10^ Mq, and mark the upper and lower limits of the measured Mo.e for 
the MW dSphs. 

centre of mas|3; 2) to assume that the inner de nsity profiles are fit 
by Einasto profiles. Following IlI et aL I l l2009l) . we compute Mo.e 
by employing the second method for satellites which show signa- 
tures of tidal perturbatiorQ and use the direct summation method 
for the other satellites. 

Our model predicts that the faintest satellites are the most dark 
matter dominated, in agreement with observational measurements. 
It is encouraging that the data points and the model satellites follow 
a similar trend, and that the spread in the mass-to-light ratio is also 
comparable. These resu lts can translate into the Mq.q-Lv plot pre- 
sented and discussed in iLi et al.l dioO^ and are consistent with the 
claim that the MW satellites are embedded in dark haloes whose 
mass i s ~ 10^ Mq within the optical extent as proposed bv lMateol 
( Il998l) . despite th e fact that their lum inosities span nearly five or- 
ders of magnitude jStrigari et alj2008h . The same conclusions hold 
when we use the luminosities predicted by the alternative feedback 
scheme, and the only difference is that there are more faint satellites 
(Mv > -5) which have 3, 000 < Mos/W < 20, 000 [M/L]q. 

Datasets which cover large radius as well as internal proper 
motions for stars in these systems are crucial for drawing further 
conclusions on the total mass content of the MW satellites. From 
our modelling of the baryonic physics, we expect a minimum dark 
matter halo mass of galaxies before accretion onto a MW-like host 
(equivalent K-ir = 16.7 km s~^), introduced by the atomic hydro- 
gen cooling limit (see also the discussion in Section |43t . 



4 DISCUSSION AND IMPLICATIONS 
4.1 Number of satellites 

Spectra of distant quasars s uggest that reionization was completed 
by 2 = 6 jFan et al1l2002[) . yet the exact duration and processes 

^ determined with the 10 per cent most bound particles in each associated 
subhalo 

* The two model satellites with Af* > A/dm at z = are excluded from 
this analysis. 



by which the Universe was reionized are not well understood. Our 
choice of 2:roio = 15 is consistent with the current observational 
constraints of Zrdo = 11-0 ± 1.4 given by the WMAP 5-year 
data, which a re obtained assuming the reionization occurred in- 
stantaneously jHinshaw et al.l2009l) . In our model, the suppression 
of cooling in haloes with Vvir < 16.7 km is the most crucial 
mechanism which alleviates the discrepancy in the number of sur- 
viving dark and luminous satellites. The number of luminous satel- 
lites is 1774 in a model that does not take reionization into account 
and in which haloes with Tvir lower than 10* K, are able to cool 
as much gas as a 10* K halo with the same metallicity. If cooling 
is forbidden in these haloes, A^'sat reduces to 121. In a model that 
includes only reionization (no suppression of the cooling in small 
haloes) A^sat ~ 286 if Zroio ~ 8, and Nsat ~ 73 when Zraio ~ 15. 
Forbidding cooling in small haloes, the number of luminous satel- 
lites varies from 88 to 51 when changing Zroio from 8 to 15. The 
dependence of A'^sat on Zrcio in our mo del appears slightly stro nger 
than in other recen t studies, e.g. by iKravtsov et al.l ( |2004|) and 
iMaccio etalj | |2009|) . Both Afsat(2reio = 15) and Afsat(zreio = 8) 
from our modelling are, however, consistent with current observa- 
tional constraints, especially when the possible anisotropic distri- 
bution of the satellites is taken into account l^llerud et al. 200^). 
The dependence on Zroio found in our model appears to be consis- 
tent with recent results bv lBusha et al] l l2009h . 

Recent hydrodyna mical simulation s of reionization have sug- 
gested that the work by lOnedinl 1 2000l) might have overestimated 
the va lue of the 'filtering mass' jHoeft et^l2006l : IOkamoto et al.l 
l2008h . We acknowledge that this could imply that our model may 
underestimate the number of surviving satellites, especially at low 
or intermediate luminosities. iMaccio et al.l ( l2009h show, however, 
that effect on the luminosity function of lowering the value of filter- 
ing mass can be compensated with an earlier epoch of reionization. 
We remind the reader that our study could not address the scat- 
ter in the luminosity fu nction expected as a result of halo-to-halo 
variation. IMaccio et alj show that the number of satellites changes 
by about a factor of two in each luminosity bin owing to the vari- 
ations of merging histories of host galaxies in the mass range of 
0.8 — 1.2 X 10^^ h"^ Mq. Their study also suggests that more mas- 
sive central galaxies (with M = 2.63 x 10^^h~^ Mq) would host 
a larger number of satellites. These studies can be addressed us- 
ing new high-resolution simulations of MW-like haloes covering 
a range of total masses and merging histories (e.g. the Aquarius 
project. ISpringel et al.l2008h . 

Our model does not take into account the loss of stars due 
to tidal stripping. We have, however, analysed the evolution of the 
bound dark matter mass (A/dm) and A/o.e for the surviving model 
satellites, and found that only ^ 10 of them show significant evolu- 
tion in the last 2 Gyrs. Since baryons are much more concentrated 
than dark matter, we do not expect this process to affect signifi- 
cantly the luminosity of our satellites. 

Another process that we have not included, and which is po- 
tentially important for the satellite luminosities is ram pressure 
stripping of the cold gas. This is likely to be a stronger effect in 
our fiducial satellite-model than in the ejection model, where SN 
feedback itself is strong enough and thus reduces the cold gas con- 
tent of satellite galaxies. The impact of ram-pressue on the most 
luminous satellites is also expected to be weak as these are typ- 
ically associated with the most massive subhaloes, that were ac- 
creted most recently. A more careful modelling that includes how 
baryons are aff ected by the interac tions with the MW-like galaxy is 
clearly needed jMaver et alj2006l) . 



4.2 Satellites associated witti small dark matter haloes and 
ultra-faint satellites 

In our study, we have excluded satellite (Type 2) galaxies whose 
associated dark matter subhalo was reduced below the resolution 
limit of the simulation (i.e. < 2 x 10*^ M0). In our fiducial satellite- 
model, there are 72 Type 2 galaxies, that cover a wide range of 
present day luminosities (similar to that covered by satellites asso- 
ciated with a distinct dark matter substructure but extended to the 
fainter end). By tracing their position using the most bound parti- 
cles of their associate subhaloes at the last time they were identified, 
we find that the median distance of Type 2 galaxies from the host is 
~ 20 kpc at z = 0. 90 per cent of them are within 120 kpc with a 
maximum distance of ~ 170 kpc, and they follow a much more cen- 
trally concentrated distribution compared to that of Type 1 galaxies 
(see also iGao et al, 2004). Among those Type 2 galaxies within 
the iimer 20 kpc, about 40 per cent are fainter than Mv = — 5 at 
z = using the more effective SN feedback scheme. If we were 
to find the counterparts/remnants of Type 2 galaxies in the observa- 
tions, they would likely show a disturbed morphology and/or tidal 
streams and reside very close to the host galaxy. A few ultra-faint 
satellites which show irregular morphologie s might be indicating 
they are ti dally disturbed, e.g. Hercules dSph jColeman et alj2007l : 
ISand et al . 2009). However, it is not yet clear whether the irregular- 
ities seen in some ultra-faints are indeed ow ing to tidal interacti ons 
given the scarce of stars currently observed l lMartin et alj2008h . 

The luminosity function predicted by our model is in good 
agreement with the power-law-shaped 'all-sky' SDSS luminosity 
function, especially when using the more efficient SN feedback 
scheme. However, this model underestimates the number of ex- 
pected ultra-faint (fainter than Mv = —5) satellites (see Fig. |4] 
right panel). Recent theoretical studies have proposed that the ultra- 
faint satellites are those fossil satellites which for med before reion- 
ization and managed to cool some gas via H 2 1 Bovill & Ricottil 
I2OO9I : ISalvadori & Ferraral |2009| ; iMufioz et alj|2009l) . In Fig. [T4l 
we show the luminosity functions with three different models all 
using our ejection SN feedback scheme. The lightest histogram 
shows results obtained using Zreio = 8 and allowing haloes with 
Tvir < lO** K to cool as much gas as a Tvir = 10"* K halo of 
the same metallicity. The intermediate histogram adopts the same 
assumption for cooling in small haloes but assumes z^cio ~ 12. 
The thickest histogram also adopts Zido = 12 but forbids cool- 
ing in small haloes. The two models using Zraio = 12 agree well 
with each other except at the faint end (as expected). So shutting 
down cooling in small haloes mostly affects the luminosity func- 
tion fainter than Mv = —5. If some of those small haloes could 
have cooled via molecular hydrogen, they would indeed populate 
the ultra-faint satellites regime. However, it should be noted that, 
when adopting Zreio ~ 12, there is still a non negligible fraction 
(about 20 per cent) of ultra-faint satellites that are associated with 
haloes with Tvir > 10* K and whose low luminosities are driven 
by the very efficient SN feedback. 



4.3 Which dark matter substructures can form stars? 

Fig.[T]shows that the mass range 10*^ - lO'' M0 is populated with 
subhaloes both with and without stars in the fiducial model. It is 
interesting to ask why some massive subhaloes have failed to form 
stars. To address this question, we trace the evolution of the dark 
matter mass M200 for all satellites and store the maximum mass 
and epoch (t = tm) when this is reached. Fig.[T5]shows the maxi- 
mum mass against the present-day (bound) mass Mom - Luminous 
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Figure 14. Luminosity functions for the MW satellites (dashed-line and 
points) and for our model satellites using the more efficient supernova feed- 
back scheme (solid histograms). Data are the same as those in Fig. [3] 
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Figure 15. Maximum A/200 mass against present-day bound dark matter 
mass for luminous (black) and dark (grey) satellites. The symbol size for the 
luminous satellites increases with the luminosity as given by the satellite- 
model. 



satellites are shown in black and dark satellites in grey. The symbol 
size for the satellites increases with the luminosity as predicted by 
the fiducial satellite-model. It is clear that most of the subhaloes 
associated with the luminous satellites were once much more mas- 
sive, and that the most luminous satellites are on average embedded 
within the most massive subhaloes at present-day (see also Fig.|6]l. 
On the other hand, dark satellites (in this mass rang e) have present- 
day m asses similar to their peak values (see also iKravtsov et al.l 
l2004h . It is also clear that model satellites have a wide range of 
peak M200 prior to their accretions onto the central galaxy. 

We plot the maximum mass as a function of the redshift when 
this value was reached in Fig. [16] The minimum mass defined by 
the cooling via atomic hydrogen as a function of redshift is indi- 
cated by the solid curve. The Figure shows that most of the dark 
satellites were below the threshold and not able to cool gas, even 
when they reached their peak mass. In contrast, luminous satellites 
live in subhaloes which have been massive enough and managed to 
have sufficient cold gas to fuel star formation. We also notice that 
no luminous satellites achieved their maximum mass before z = Q. 
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Figure 16. Maximum M200 mass as a function of the redshift when it was 
reached for luminous (black) and dark (grey) satellites. The minimum mass 
for cooling via atomic hydrogen at each redshift is indicated by the solid 
line. Symbols are the same as in Fig. 1151 



We have also checked when and where the present-day lumi- 
nous sateUites typically formed in the simulation. We find that the 
surviving satellites were all detected in the simulation early, namely 
their Mom exceeded 2 x 10*' M0 around the epoch of reionization 
(zreio ~ 15 in the satellite-model) and mostly before the end of 
reionization (z = 11), while the dark ones have emerged through- 
out the Hubble time. The distance of the model satellites to the 
MW -like galaxy at em ergence is typically around ~ 100 kpc. 

lFontetaiT ( l2006l) suggest that surviving satellites were ac- 
creted up to 9 Gyr ago and therefore have been through a different 
chemical enrichment history compared to those that were accreted 
very early on and contributed to the stellar halo. In our model there 
are a few objects that became satellites more than 10 Gyr ago and 
have survived the tidal interactions with the MW-halo (see Fig.|6] 
and Fig.[8]l. Typically, however, their peak masses are smaller than 
those of the object s that contributed signific antly to the build up 
of the stellar halo toe Lucia & Helmill2008l) . It is possible to en- 
vision the building blocks of the stellar halo as the most massive 
satellites that existed at early times, i.e. the counterparts of the very 
faint ones, whose orbits did not decay via dynamical friction, and 
hence have generally been more sheltered from the tidal forces of 
the Galaxy. 



5 CONCLUSIONS 

We use a hybrid model of galaxy formation and evolution to study 
the satellites of the Milky Way in a cosmological context. Our 
method combines high resolution A'^-body simulations which allow 
us to trace the evolution and the dynamics of dark matter haloes 
directly, and phenomenological prescriptions to follow the evolu- 
tion of baryons. Our adopted semi-analytical recipes and values for 
the relevant parameters result in models that reproduce the prop- 
erties of galaxies on large scales as well as those of the MW. A 
few modifications were needed, however, to reproduce the obser- 
vational properties of the MW satellites. 

With the presence of a reionization background that reduces 
the baryon content of subhaloes around z = 15 and the suppression 
of cooling for haloes with Vvir < 16.7kms~^, our model can re- 
produce the total number and the luminosity function observed for 



the satellites of the MW. Our fiducial SA model also shows good 
agreement with the metallicity distribution and the metallicity- 
luminosity relation when a large fraction of newly formed metals 
in small galaxies are recycled through its hot component and with 
other properties shared by the MW satellites, e.g. the radial dis- 
tribution, luminosity-size relation and the star formation histories. 
However, our fiducial model produces an excess of satellites with 
AIv ~ —10 and log(Z,/ Z©) ~ —1, and does not predict ultra- 
faint satellites with the total luminosity below Lv ~ 10*'^ Lq. 

We have tested an alternative SN feedback recipe which is 
stronger for galaxies with Vvii ^ 90kms~^ compared to the stan- 
dard feedback recipe. With this alternative feedback recipe, our 
model predicts the same number of surviving satellites (which pop- 
ulate the same set of subhaloes as the standard feedback model). 
The alternative feedback model predicts more satellites with My < 
—5 which also follow the metallicity-luminosity relati on traced by 
the cl assical and the ultra-faint SDSS MW satellites jKirbv et"al] 
|2008|) down to [Fe/H] ~ -2.5 and Lv ~ 10^ Lq. 

Our model satellites are embedded in dark matter haloes 
with innermost masses within 600 pc between 6 x 10^ M0 and 
7 X 10^ M0, in very good agreement with the estimates for the 
classical MW dSphs derived by Strigari et al. (2003). This demon- 
strates that the existence of a common scale for the innermost mass 
is a natural outc ome of the CDM galaxy formation and evolution 
mode l (see also iLi et al.|[2009l ; iMaccio et Z||2009| ; IKodosov etal] 
I2OO9I) . Satellites are dark matter dominated even within the optical 
extent, and the mass-to-light ratio increases with decreasing lumi- 
nosity. Surviving satellites in our model are associated with ancient 
haloes which had masses of a few 10** Mq by z ~ 10 - 20, and 
acquired their maximum dark matter mass after z ^ 6. 

The brightest satellites in our model are associated with the 
most massive subhaloes, were accreted later (z < 1) and show 
extended star formation histories, with only 1 per cent of their stars 
made by the end of the reionization (lookback time ~ 13 Gyr). On 
the other hand, the fainter satellites tend to be accreted early on, are 
all dominated by stars with age > 10 Gyr, and a few of them are 
dominated by stars formed before the reionization was complete. 
In our models, the classical MW satellites are associated with dark 
matter subhaloes with a peak circu lar velocity > 10k ms~^, in 
agreement with the recent results bv lWalker et al] I2OO9I) . 

Although our model satellites are in very good agreement with 
the latest observations in terms of luminosity, metallicity and the 
innermost dark matter content, the agreement with the ultra-faint 
galaxies is less conclusive. Note that the alternative feedback recipe 
seems to give a better fit to the properties of ultra-faint satellites. 
This could imply that the ultra-faint satellites of the MW are as- 
sociated with dark matter haloes with lower Vvir compared to the 
classical ones (i.e. they are more sensitive to SN feedback). How- 
ever, until we properly model the loss of baryons due to tidal strip- 
ping and ram pressure in these small systems, and until we obtain 
more observational constraints on these objects, their nature will 
probably remain unclear. 
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